On the Lieb-Liniger model in the infinite coupling constant limit 
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We consider the one-dimensional Lieb-Liniger model (bosons interacting via 2-body delta poten- 
tials) in the infinite coupling constant limit (the so-called Tonks-Girardeau model). This model 
might be relevant as a description of atomic Bose gases confined in a one-dimensional geometry. It 
is known to have a fermionic spectrum since the A''-body wavefunctions have to vanish at coinciding 
points, and therefore be symmetrizations of fermionic Slater wavefunctions. We argue that in the 
infinite coupling constant limit the model is indistinguishable from free fermions, i.e., all physically 
accessible observables are the same as those of free fermions. Therefore, Bose-Einstein condensate 
experiments at finite energy that preserve the one-dimensional geometry cannot test any bosonic 
characteristic of such a model. 
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Statistics in one dimension and the hard-core 
boson model: The Lieb-Liniger model of interacting 
particles on the line is defined as [1] 
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Contact (^-interactions have no effect on fermionic wave- 
functions whereas they act nontrivially on bosonic ones. 
Such interactions are commonly used to model interact- 
ing atomic bosonic gases [2] and can also be relevant in 
describing other aspects of atomic quantum fluids such 
as, for example, supersolidity Q. 

The Lieb-Liniger model is solvable by Bethe Ansatz 
techniques for any finite coupling c. In the c — s- cx) hard 
core limit, the spectrum of the model becomes fermionic 
since the eigenstates have to vanish at coincidence points. 
On the other hand the statistics of the model, conven- 
tionally defined through the symmetry of the wavefunc- 
tion, remains bosonic. The energy eigenfunctions V'n.G 
are symmetrizations of the fermionic ones ipn,F, 

i<j 

(2) 

and the energy eigenvalues are identical to the fermionic 
ones. This particular Bqse-Fermi mapping has been first 
noticed by Girardeau |4)]. The question we address in 
this note is the exact fermionic nature of the Lieb-Liniger 
model in this singular limit. 

In one dimension, the configuration space for N par- 
ticles decomposes into A^! distinct sectors, according to 
the ordering of the particle coordinates. If the particles 
are identical and indistinguishable, these sectors are all 
physically equivalent. In principle we could keep only one 
sector, the other ones representing "gauge" (unphysical) 
copies. 

Consider an A'^-body quantum Hamiltonian and an ini- 
tial wavefunction with support inside a particular sector 



of the configuration space, say, xi < X2 < ■■■ < xn. For a 
Hamiltonian without singular interactions, the quantum 
time evolution of the wavefunction is such that differ- 
ent sectors end up communicating with each other; that 
is, the wavefunction will eventually spread over the re- 
maining sectors. The different sectors interfere and it is 
necessary to determine the value of the wavefunction in 
all of them to fully fix the dynamics. 

However, when the interaction becomes sufficiently sin- 
gular at coincidence points, tunneling between different 
sectors might become forbidden, implying that the wave- 
function time evolution is restricted to the initial sector. 
In such a case, wavefunctions from different sectors never 
interfere and live independent lives. The sectors have be- 
come effectively superselected and particle statistics have 
become irrelevant: we may restrict the wavefunction in 
one sector, or continue it in a symmetric or antisymmet- 
ric way in the other sectors without affecting the physics. 

This situation is precisely encountered in the 
Calogero model where the singular nature of the 
g{g — l)/{xi — XjY Calogero interaction prohibits quan- 
tum tunneling between different sectors which therefore 
do not communicate. The inverse square potential is 
quantum mechanically impenetrable in such a way that 
if the particles are in a given sector among the A^! possi- 
ble, they will stay in this sector for ever. This property 
of the Calogero interactions is an important component 
in understanding the excusion statistics pertaining to the 
model [Hi. 

In the hard core Lieb-Liniger model the situation is 
analogous: because of the singular nature of the 5 in- 
teractions in the infinite coupling constant limit particle 
penetration is suppressed. Indeed, the tunneling ampli- 
tude for a particle of momentum p scattering off a delta- 
potential of strength c is 



ip 



ip - 



(3) 



2 



and tunneling for momenta lower than c is suppressed. 
In the limit c — > oo all tunneling is suppressed. 

In the reference sector xi < X2 < ■■■ < xn the 
hard core Lieb-Liniger wave functions ^l>G coincide with 
fermionic Slater wavefunctions ipp- Since the original 
statistics of the model is assumed to be bosonic, one 
may continue the wavefunctions to other sectors of the 
configuration space in a symmetric way as in ([2]) • But 
this continuation is physically irrelevant since, as already 
stated, tunneling to the other sectors is impossible. What 
matters for physics is the wavefunction in the reference 
sector. As this is the same as that of free fermions, the 
system is entirely fermionic. It follows that the model, 
which decribes free bosons when c = 0, ends up describ- 
ing free fermions in the c — s- oo limit. 

For exchange or fractional statistics models |6j , such as 
the two-dimensional anyon model or the one-dimensional 
Calogero model, statistical interactions are parametrized 
by a dimensionless statistical coupling constant -in the 
anyon model it is the fractional part of the flux carried by 
the particles in unit of the flux quantum; in the Calogero 
model it is g. This absence of a scale is natural, as one 
does not expect statistical considerations to arise from a 
dimensionful interaction, which would introduce an ad- 
ditional physical scale. In the Lieb-Liniger model the 
coupling constant c has dimensions of inverse length, so 
it is not a statistical model per se except when trivially 
c = 0, i.e. pure bosons, and when less trivially c = oo, 
i.e. pure fermions. At both points the dimensionful scale 
c indeed disappears from the model. In between one can 
still view the Lieb-Liniger model as some sort of statisti- 
cal model with an energy-dependent statistical parame- 
ter, which for energies much higher than is essentially 
bosonic, while for energies much lower than is essen- 
tially fermionic. 

An alternative point of view for the statistics of par- 
ticles in one spatial dimension is to restrict the wave- 
function from the outset in one sector (a kind of 'gauge 
fixing'). Since each sector has a boundary (where any 
two particle coordinates coincide) the issue of bound- 
ary conditions on the wavefunction becomes relevant 0] . 
Neumann (9„i/'6 = 0) or Dirichlet {iph — 0) are obvious 
choices that preserve hermiticity and are equivalent to 
symmetric (bosonic) or antisymmetric (fermionic) con- 
tinuations of the wavefunction over the remaining sec- 
tors, respectively. Linear combinations of the above of 
the sort dnipb + cij^b = are also possible, and correspond 
to the bosonic Lieb-Liniger model of strength c. Once 
again, we recover the case c = cxd as corresponding to 
fermionic statistics. 

Physical observables and momentum: One cru- 
cial aspect that makes the Tonks- Girardeau gas appar- 
ently distinct from free fermions is its momentum dis- 
tribution. As the wavefunctions of the two systems dif- 
fer only by their phase (sign) in each sector, coordinate 
observables have the same expectation value, but the 



momentum wavefunctions (Fourier transforms) are gen- 
uinely different Q. How can this be reconciled with the 
statement that the two models are physically equivalent? 

To answer this we must address the question of 
what are physically accessible observables in the Tonks- 
Girardeau system. Concentrating on the momentum, any 
local interaction of the system with a measuring device 
would involve a fully symmetric polynomial of finite de- 
gree in the particle momenta pi, . . . ,pN. 

To reduce the problem to its bare bones, we examine 
the simplest stituation of two particles. The center of 
mass momentum pi +p2 is an observable, as well as all its 
powers. The relative momentum p — [pi — P2)/2, on the 
other hand, is not an observable since it is not invariant 
under particle permutations, but its even powers p^" are. 

It is now easy to establish that the center of mass 
momentum has the same matrix elements in both mod- 
els. Indeed, since the particle interaction in the Tonks- 
Girardeau model is translation invariant, the center of 
mass motion is free. We may factor out the center of 
mass motion and consider only the relative part, with 
wavefunction ip{^)j ^ = xi — X2- For the relative mo- 
mentum observables p^" we have 

<p2">= J da;7/>*(-i9^)2"V = J dx\d^^\^ (4) 

where we integrated by parts n times and assumed the 
usual regular behavior oi tp ai x —^ ±oo. Since ijjGix) = 
sgn{x)ijjF{x) and both wavefunctions vanish linearly at 
X = we have 

d^ijG = sgn(x)d,i;F , d^^^a - S^"-^\x) , n>2 (5) 

That is, ipG has a discontinuity in the first derivative 
and thus develops as the leading singular term a delta 
function of order n — 1 at a; = in the 7i-th derivative. 
^From (|4]) and (O we conclude that 

<P^>G^<P^>F , < p^" >G= OO , 71 > 2 (6) 

Therefore, all observables that would be different be- 
tween the two models have infinite expectation value in 
the Tonks-Girardeau model and thus are not physically 
measurable, requiring an infinite energy to do so. The 
physical reason behind this effect is that such operators, 
added as perturbations to the system's Hamiltonian due 
to the coupling to the measuring device, have the ef- 
fect of "pushing" particles through each other and the 
hard-core potential between particles imposes an infinite 
energy penalty to this process. 

In a similar way we can show that mixed local op- 
erators of the form 2;™p^"~™ (or their hermitian syni- 
metrization) follow a similar pattern: 

< >G = < >^ ^ n - m < 1 

< X V""" >G = oo , n - m > 2 (7) 
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That is, operators that are not the same as in the fermion 
system have infinite expectation value and are, thus, in- 
accessible. 

There are, of course, other p-dependent operators that 
are different between the two models and do not have 
infinite expectation value, such as, e.g., cos(ap) with a 
a constant. These operators, however, are nonlocal and 
have the property to "teleport" particles over each other, 
thus avoiding the interaction region and the correspond- 
ing energy cost. (The two-particle permutation opera- 
tor itself is another example.) Such observables are not 
physically realizable. We can argue that operators mix- 
ing the superselection sectors manipulate gauge copies of 
the configuration space and are, in general, unphysical. 

Time evolution and momentum operator: From 
the discussion of the preceding section it becomes clear 
that operators such as p"* become unphysical for the 
Tonks-Girardeau model and therefore the momentum op- 
erator itself is unsuitable. The problem is similar to the 
quantum mechanical description of a particle on the half- 
line. It is known that the standard operators x and p are 
unsuitable, and a minimal set of physical operators is 
x^, and xp + px, closing into the SL{2) algebra. The 
Casimir of this algebra fixes the physical Hilbert space 
and implies that p may not be represented as —idx on 
this space. 

To identify an appropriate momentum operator we pre- 
fer here to follow a more physical approach and go back to 
basics, remembering the justification for defining —idx in 
the standard Schrodinger picture as the momentum op- 
erator or, equivalently, the Fourier transform of the co- 
ordinate wavefunction as the momentum wavefunction. 
We shall work again with the simplest case of two parti- 
cles and factor out the center of mass coordinate, leaving 
us with the relative coordinate x. 

Consider a free system with initial wavefunction in the 
vicinity of a: = and an overall uncertainty in the po- 
sition X of order Aa;. We let the system evolve for a 
long time T and measure the position x. Typically x 
will have some large value L, as the particles have moved 
away from each other. If we measure x to be in the 
range [L — AL, L + AL] we know that it moved a dis- 
tance D = X f — Xi = L± AL ± Ax in time T. Its speed, 
and correspondingly the momentum (assuming unit re- 
duced mass), is estimated then as 



tion operator 



L AL Ax 
n = V = — ± ± 

T' T T 



(8) 



By taking T large enough, and scaling L and AL with 
T, we can make the uncertainty of p due to the ini- 
tial Aa; negligible. The probability of the momentum, 
then, being within a range Ap around the value p is the 
probability of measuring the coordinate a; within a range 
AL = TAp around the value L = Tp, for T — )■ oo. 

The large-T evolution of the wavefunction can be cal- 
culated from the initial wavefunction and the free evolu- 



U{T) = e 



(9) 



In terms of the Fourier transform of the wavefunction 
(j){k) we have 



^(a;,T) 



1 



(10) 



Putting X = Tp and performing a standard saddle-point 
expansion for large T we obtain 



^(rp,r) = -^0(p) (T^oo) 



(11) 



Since the probability for measuring x in a range TAp 
around Tp is TAp\'ij:{Tp, T)p, we find the probability for 
measuring the momentum in the range Ap around p to 
be 



P{p,P + Ap) ^ Ap\(l>{p)[' 



(12) 



which justifies the role of the Fourier transform ^(p) as 
the momentum wavefunction. 

We can now use this procedure to define the momen- 
tum of particles in the Tonks-Girardeau system: let the 
particles fly away from each other. The distribution of 
their positions after they spread enough will be taken 
to represent their kinematical velocity, or momentum, 
distribution. Although this evolution is not really free, 
interactions occur only at coincidence points and are an 
integral part of the definition of the system, just like the 
restriction to a; > for the particle on the half-line. 

It should be clear that this distribution will be iden- 
tical to the one of freely moving fermions. To see this, 
use the fact that the spectra of the two models are iden- 
tical and that the energy eigenfunctions are related as in 
([2]). Then, expanding the initial wavefunction in terms 
of energy eigenstates, we have 

iIj^{xi,...,Xn) = ^Cn1pn,G{xi,--,XN) 

n 

= Ck TT sgn{xi - Xj)ipn,F{xi, Xn) 

n i<j 

= Wsgn{x^ - Xj) Cnll)n,F{xi, ■■■,Xn) 
i<j n 

= Y\_^9'n-ixt - Xj)TpFixi,...,XN) (13) 

In other words, the hard-core wavefunction and the cor- 
responding fermionic wavefunction have the same energy 
eigenvalue expansion coefficients. Thus, at any time t, 

ip^ixi, ...,XN,t) = ^c„e"*^"*'0„,G(a;i, ...,a:Ar) 

n 

= Y\_sgn{xi- Xj)tpF{xi,...,XN,t) 

(14) 
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So the wavefunctions of the hard-core boson model and 
the free fermion model remain at all times related by a 
simple sign change over sectors. Therefore the particle 
distribution after a large time T, and so the induced mo- 
mentum distribution, are that of free fermions. 

For general interacting systems, one can recover the 
standard momentum distribution by switching off the in- 
teractions and then letting the particles fly away. This 
would reproduce the (bosonic) Fourier transform as the 
momentum distribution. The arguments of the preceding 
sections assume the permanent presence of the infinite- 
strength interactions. Our point is that any real fly-away 
experiment preserving the one- dimensional geometry and 
interactions would gather data consistent with an initial 
free fermion momentum distribution. 

The above can also be summarized by saying that a 
complete set of physical observablcs is the set of all time- 
dependent amplitudes of the form 

|< xi,. .. ,XN;t\x[, . . ■,x'j^;t' >f = \G{x[ ~ Xi\t! - t)f 

From the present analysis it is clear that these amplitudes 
for the Tonks- Girardeau model are identical to those of 
free fermions. 

Conclusions: We have argued that the Tonks- 
Girardcau model is physically indistinguishable from free 
fermions, on the basis that any experiment aiming to 
distinguish the two models while staying within the one- 
dimensional setting would require an infinite energy or 
the action of unphysical operators. Inferring particle mo- 
menta from the fly-away spatial distribution of Tonks- 
Girardeau particles would return the distribution of free 
fermions. 

The above considerations are, in fact, quite relevant 
to the experimental setup decribed in [91, llO|] : allow the 
Tonks-Girardeau atoms, originally constrained in a trap, 
to fly away ballistically while they are still in a Id ge- 
ometry and use the asymptotic particle distribution to 
deduce the momentum distribution. As argued above, 
such an experiment would fail to detect any bosonic fea- 
tures of the system. 

In experimental realizations of the model the interac- 
tions are not entirely hard-core; that is, the strength of 
the potential will not be infinite. For a finite c the tunnel- 
ing between sectors is not suppressed for momenta p > c 
and the energy cost of measurements that would dist- 
inghuish the model from fermions is not infinite. In effect, 
the model describes free fermions for momenta p <^ c and 
free bosons for momenta p ^ c, with a crossover region 
in between. Interpretation of experimental data, there- 
fore, must take into account the finite strength of the 
core repulsion relative to the momentum scale. 



The arguments presented in the previous analysis can 
be repeated and verified in the discrete lattice version of 
the model, that is, the Bose-Hubbard model [HI with 
infinite on-site repulsion energy and finite hopping pa- 
rameter. In this model at most one boson is allowed on 
a given site, so again the particles are fermion-like. The 
bosonic and fermionic Hilbert spaces are in exact one-to- 
one correspondence via a Jordan- Wigner transformation 
that maps each bosonic operator into a fermionic one, 
and vice versa. The total momentum operator is com- 
mon in the two descriptions, but the relative momentum 
operators differ: the bosonic one may create states of in- 
finite energy (since it may move two bosons on top of 
each other), while the fermionic one does not. Again, we 
recover the same basic conclusions as in the continuous 
model. 
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